We derive accurate parameters related to the CMD, 
structure and dynamical state of M52 and NGC3960, 
whose fields are affected by differential reddening. 
Previous works estimated their ages in the ranges 35 — 
135 Myr and 0.5 — 1.0 Gyr, respectively. J, H and K s 
2MASS photometry with errors < 0.2 mag is used to build 
CMDs, radial density profiles, luminosity and mass func- 
tions, and correct for differential reddening. Field-star de- 
contamination is applied to uncover the cluster's intrinsic 
CMD morphology, and colour-magnitude filters are used 
to isolate stars with high probability of being cluster mem- 
bers. The differential-reddening corrected radial density 
profile of M 52 follows King's law with core and limiting 
\Q radii of R COTe = 0.91 ± 0.14 pc and R lim = 8.0 ± 1.0 pc. 

■ NGC 3960 presents an excess of the stellar density over 
, King's profile (R CO ic = 0.62 ± 0.11 pc and i?i; m = 6.0 ± 

^ 0.8 pc) at the center. The tidal radii of M 52 and NGC 3960 
are i?tidai = 13.1 ± 2.2 pc and i?tidai = 10.7 ± 3.7 pc. 
£3 ' Cluster ages of M 52 and NGC 3960 derived with Padova 
■^-s , isochrones are constrained to 60±10Myr and 1.1 ±0.1 Gyr. 

■ In M52 the core MF (xcore = 0.89 ± 0.12) is flatter than 
the halo's (xhaio = 1-65 ± 0.12). In NGC 3960 they are 

. Xcorc = _o.74 ± 0.35 and Xhaio = 1-26 ± 0.26. The mass 
^ locked up in MS/evolved stars in M 52 is ~ 1200 M , and 
L^l ' the total mass (extrapolated to 0.08 M®) is ~ 3800 Mq. 
q The total mass in NGC 3960 is ~ 1300 M Q . Compared to 
OO ' open clusters in different dynamical states studied with 
, similar methods, the core and overall parameters of M 52 
^§ are consistent with an open cluster more massive than 
1 000 Mq and ~ 60 Myr old, with some mass segregation 
in the inner region. The core of NGC 3960 is in an ad- 
vanced dynamical state with strong mass segregation in 



o 



Q ■ the core/halo region, while the somewhat flat overall MF 
,H '. (y w 1.07) suggests low-mass star evaporation. The excess 
& ' stellar density in the core may suggest post-core collapse. 
. . | The dynamical evolution of NGC 3960 may have been ac- 
celerated by the tidal Galactic field, since it lies ks 0.5 kpc 
inside the Solar circle. 
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Abstract. 

1. Introduction 

Open clusters (OCs) can be used as test-beds of molec- 
ular cloud fragmentation, star formation, and stellar and 
dynamical evolution models. They are excellent probes of 
the Galactic disc structure (Janes & Phelps 119941 Friel 
1995 Bonatto et al. 2006a). This follows from the relative 
simplicity in age and distance estimates, at least for bright 
clusters. Accurate OC parameters such as core and limit- 
ing radii, age, mass, density, and relaxation times are es- 
sential to constrain theoretical models. However, because 
most of the OCs lie close to the disk, and thus are more 
affected both by reddening and field-star contamination 
(Bonatto et al. '20 06al and references therein), this kind 
of analysis becomes restricted to the more populous OCs 
and/or those located a few kpc from the Sun. 

According to the WEBDA 1 OC database (Mermilliod 
1996} the current number of OCs with known parameters 
(such as coordinates, age, distance from the Sun, and red- 
dening) amounts to ~ 700. The catalogue of Dias et al. 
(2002) includes 1756 optically visible OCs and candidates. 
More recently Bonatto et al. (2006a) simulated observa- 
tional completeness in different directions throughout the 
Galaxy and estimated that ~ 50% of Trumpler types I 
to III OCs remain undetected in the solar neighbourhood 
(distance from the Sun c? Q < 1.3 kpc). A large fraction of 
the faint and/or poorly populated OCs end up drowned 
in the stellar background, particularly in directions inter- 
cepting the bulge. 

Several decades of studies have provided fundamen- 
tal parameters for a considerable fraction of the cata- 
logued OCs (e.g. Lynga 115521 WEBD A; Dias et al. 120021 
Kharchenko et al. 120051 Bonatto et al. l2006a|l . Some clus- 
ters have been studied twice or more and a disagree- 
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ment of a factor of ~ 2 in parameters is not uncom- 
mon. Discrepancies arise mostly from different analytical 
methods dealing with optical or IR photometry, Colour- 
Magnitude Diagrams (CMDs), isochrone sets, photomet- 
ric systems, etc. 

Recently we started a program to derive fundamental 
parameters of OCs with near-IR photometry and a series 
of analytical procedures developed to minimise field-star 
contamination of CMDs. They have proven to be effective 
for parameter determination for a wide class of OCs. For 
spatial and photometric uniformity we work with J, H 
and K s 2MASS 2 photometry. The 2MASS Point Source 
Catalogue (PSC) is uniform reaching relatively faint mag- 
nitudes covering nearly all the sky, allowing a proper back- 
ground definition for clusters with large angular sizes (e.g. 
Bonatto, Bica & Santos Jr. 120051 Bonatto, Bica & Pavani 
12004) . 

The use of field-star decontamination and colour- 
magnitude filters have produced more robust parame- 
ters (e.g. Bonatto & Bica 120051 Bonatto et al. I2006bjl . 
In particular, field-star decontamination constrains the 
age more, especially for low-latitude OCs (Bonatto et al. 
2006a). These procedures have proven useful in the anal- 
ysis also of faint and/or distant OCs (Bica, Bonatto & 
Dutra l2003l Bonatto, Bica & Dutra 12001 Bica & Bonatto 

EES- 

Another potential source of observational uncertainties 
in OC parameters is the presence of differential reddening. 
Spatial variation of reddening may affect the age determi- 
nation by isochrone fitting and, consequently, the distance 
from the Sun. This in turn could affect the absolute values 
of core and limiting radii, and cluster mass. 

Our goals in this paper are twofold: (i) add to the 
previous techniques a procedure to take into account the 

2 The Two Micron All Sky Survey, All Sky 
data release (Skrutskie et al. I1997II . available at 
http: //www. ipac. caltech. edu/2mass /releases /allsky/ 
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Fig. 1. 20' x 20' XDSS R images of M52 (left panel) and NGC 3960 (right panel). 



differential reddening, and (ii) derive more constrained pa- 
rameters for two OCs of different ages. As targets we se- 
lected the relatively bright OCs M 52 and NGC 3960 that 
in the literature have uncertain age determinations. The 
age of M 52 has been estimated in the range 35 — 135 Myr 
(Sect. EU, wnile that of NGC 3960 varies from 0.5 to 
1.0 Gyr fSect. l2~2"|l . Previous works have shown evidence 
of differential reddening in the fields of both clusters (e.g. 
Pandey et al. 120011 Prisinzano et al. 12004(1 . Besides red- 
dening, age and distance from the Sun, in this paper we 
derive the core and limiting radii, and the stellar mass 
stored in the core and halo subsystems. The latter param- 
eters for NGC 3960 are derived for the first time. 

This paper is organized as follows. In Sect.Elwe review 
previous findings on M52 and NGC 3960. In Sect. |3]we 
present the 2MASS data, subtract the field-star contami- 
nation, estimate the differential reddening, derive funda- 
mental cluster parameters, and analyse the cluster density 
structure. In Sect.0]we derive luminosity and mass func- 
tions (LFs and MFs), and compute stellar content prop- 
erties. In Sect. |S]we discuss the dynamical states of both 
OCs. In Sect. HJ we use diagnostic diagrams to analyse the 
structure and dynamical state of M 52 and NGC 3960 as 
compared to other OCs. Concluding remarks are given in 
Sect.0 

2. The target open clusters M 52 and NGC 3960 

2.1. M52 

M 52 (NGC 7654) is a prominent OC located in the di- 
rection of Cassiopea, in an interarm region (Fenkart & 
Schroder 1985 ). Its field is affected by a significant amount 
of variable reddening (e.g. Danford & Thomas I1981j) and 



field star contamination (Sect. I3.2JI . consistent with its 
low galactic latitude (Table The supergiant (SG) star 
BD + 60°2532 of spectral type F7Ib with V = 8.22 and 
B — V — 1.16 that is projected « 3.1' from the cluster cen- 
ter (Sect.EJ) is a probable cluster member (Schmidt[1984). 
Other bright stars are projected in the field of M 52. These 
features appear in the 20' x 20' DSS 3 R image (left panel 
of Fig.^). Below we summarise previous works on M52. 

— One of the earliest studies of M 52 is that of 
Pesch H1960fl that using UBV photolectric photometry 
found a non-uniform reddening varying in the range 
E{B - V) = 0.51 - 0.81, a distance from the Sun 
e?0 = 1.66 kpc and suggested an age similar to that 
of the Pleiades (~ 135 Myr, according to WEBDA). 

— Lindoff ifTUrJSjl derived an age of w 35 Myr for M52. 

— Schmidt (1977) using uvby-/5 photometry found 
E(B -V) = 0.45 - .53. 

— Danford & Thomas (1981} using uvby-/3 photometry 
found <i© = 1.47 kpc, a mean E(B — V) = 0.57 and 
non- uniform reddening with E(B — V) in the northern 
half of the cluster field in excess of 0.11 with respect 
to the southern half. They found an age intermediate 
between those of the a Persei cluster (~ 70 Myr) and 
the Pleiades. 

— Bruch & Sanders l|1983[) applied an empirical calibra- 
tion to convert relative to absolute mass and estimated 
a cluster mass of « 518 M©. 

— Lynga (1987) classified M52 as Trumpler type II2r, 
and measured an angular diameter of D = 12'. 



3 Extracted from the Canadian Astronomy Data Centre 
(CADC), at http://cadcwww.dao.nrc.ca/ 
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Table 1. General data on M52 and NGC 3960 



WEBDA Present work 

Cluster a(2000) 5(2000) Age E(B - V) d a(2000) 5(2000) t b Age E(B - V) Q d Q d GC d GC 

(hms) (<"") (Myr) (kpc) (hms) ("'") (°) (°) (Myr) (kpc) (kpc) (kpc) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) 

M52 23:24:31 +61:35:29 58 0.65 1.42 23:24:42 +61:35:42 112.81 +0.44 60+ 10 0.58g;|B 1.4 8.7 7.9 

NGC3960 11:50:33 -55:40:24 664 0.30 2.26 11:50:38.9 -55:40:48 294.38+6.18 1100 + 100 0.13g:o3 1.7 7.5 6.7 
Table Notes. E(B — V) c is the colour excess of the cluster's central region (Sect. 15. 5L while the bounds refer to the range of 
reddening throughout the cluster field. Uncertainty in E(B — V) Q is of the order of ±0.05. Col. 14: g!gc calculated using 
Ro = 8.0 kpc (Reid 1993| as the distance of the Sun to the Galactic center. Col. 15: dec using Ro = 7.2 kpc (Bica et al. 
2006). Uncertainties in d0 and dec are of the order of ±0.2 kpc. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



(J-H) (J-H) 



Fig. 2. J x (J-H) CMDs of the central 10' of M 52 (left panel) and NGC 3960 (right panel). Solid lines: fiducial lines 
used to derive the differential reddening (Sect. 13. 1[) . Red stars belong mostly to the field (Sect. 13.2(1 . Error bars are 
shown in the right panels of Fig. 0] 



— Kaltcheva ( 1990J using uvby photometry found an 
age of « 96 Myr and mean cluster reddening of 
E(B - V) = 0.57. 

— Battinelli, Brandimarti & Capuzzo-Dolcetta l| 19940 us- 
ing integrated UBV photometry found an angular di- 
ameter of 12', age « 35 Myr, E(B - V) = 0.57, 202 
member stars and mass m w 440 Mq . 

— Lotkin l|1994[l estimated a main sequence (MS) turn-off 
age of ~ 50 Myr and Jq — 1.49 kpc. 



- Pandey et al. (2001) using wide field (40' x 40') 
CCD UBVIc photometry found variable reddening 
(E(B-V) = 0.46 - 0.80), d Q = 1.38 ± 0.07kpc, 
and a spread in age from 30 to 100 Myr. For stars 
with mass in the range 1.5 — 4.0 Mq they found MF 
(<j>(m) — ^ oc m- {1+ ^) slopes \ = 1-07 ± 0.08 (in- 
ner region), \ = 1.28 + 0.20 (intermediate region), and 
X = 1-40 ± 0.07 (whole cluster). They concluded that 
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the presence of mass segregation is consistent with the 
cluster's age range. 

— Nilakshi et al. fl2002) using photometry from 
Palomar Observatory Sky Survey I plates found a 
Galactocentric distance dec = 9.2 kpc, core radius 
-Rcore = 1-59 ± 0.08 pc and core stellar density p c = 
22 A ± 1.7 stars pc -2 , halo radius Rh = 4.48 pc and 
halo stellar density pn — 8.4 ± 0.3 stars pc~ 2 . 

— Kharchenko et al. (|2005) derived an age of « 69 Myr, 
E(B — V)= 0.65, d Q = 1.42 kpc, i? CO rc = 4.8' (» 2pc) 
and cluster radius of 10.2' (« 4.2 pc). 

In WEBDA the central coordinates of M 52 are (J2000) 
a = 23 h 2A rn 31 s , and 6 = +61°35'29". However, the cor- 
responding radial density profile built with stars having 
colours consistent with the cluster's (Sect. 13.4(1 presented 
a dip at R = 0'. Consequently, we searched for a new 
center by examining histograms for the number of stars 
in 0.5' bins of right ascension and declination. The co- 
ordinates that maximize the central density of stars are 
(J2000) a = 23' l 24 m 42 s , and 6 = +61°35'42". In what 
follows we refer to these optimized coordinates as the cen- 
ter of M52. Additional cluster parameters from WEBDA 
are given in Tabled 

2.2. NGC 3960 

NGC 3960 is a moderately-rich OC projected in the 4th 
quadrant (Table QJ. The 20' x 20' DSS R image is in Fig. [2 
(right panel). Previous works on NGC 3960 are listed be- 
low. 

— Using BV and DDO photometry Janes 1(1981(1 esti- 
mated the age of NGC 3960 as similar to that of the 
Hyades (0.5 - 1 Gyr), a mean reddening of E(B — V) = 
0.29 ± 0.02 and a metallicity [Fe/H] = -0.30 ± 0.06. 

— Carraro, Ng & Portinari |(Tj2iB| derived an age of 
0.6 Gyr and [Fe /H] = -0.34. 

— Prisinzano et al. ( 2004) with BVI photometry of a 34' x 
33' field centered on NGC 3960 found an age in the 
range 0.9 - 1.4 Gyr, variable reddening E(V — I) = 
0.21 - 0.78 (E(B - V) = 0.17- 0.62) with E(V - 1) = 
0.36 (E(B - V) = 0.29) at the cluster center, and an 
MF slope x = 1-81 ± 0.84. 

— Bragaglia et al. (20061 used UBVI CCD photometry to 
derive an age in the range 0.6 - 0.9 Gyr, d© = 2.1 kpc, 
[Fe/H] = -0.12 ± 0.05, E(B - V) = 0.29 ± 0.02 with 
a differential reddening AE(B - V) = ±0.05. They 
detected mass segregation and dynamical evaporation. 

— Similarly to M 52 we had to recalculate the cluster's 
central coordinates with respect to those in WEBDA 
(Tabled). 

3. The 2MASS photometry 

VizieR 4 was used to extract J, H and K s 2MASS pho- 
tometry in a circular area with radius R = 40' centered 

4 http://vizier.u-strasbg.fr/viz-bin/VizieR?-source=II/246 



on the optimized coordinates of both clusters (Tabled- 
As photometric quality constraints the extraction was re- 
stricted to stars (i) brighter than the 99.9% Point Source 
Catalogue Completeness Limit 5 , J = 15.8, H = 15.1 and 
K s = 14.3, respectively, and (ii) with errors in J, H and 
K s smaller than 0.2 mag. For reddening transformations 
we use the relations Aj/Ay = 0.276, Ah /Ay = 0.176, 
A Ks /A v = 0.118 , and Aj = 2.76 x E(J - H) (Dutra, 
Santiago & Bica 12002(1 . assuming a constant total-to- 
selective absorption ratio Ry =3.1. 

The J X (J-H) CMD of the central 10' of M52 is 
in Fig. [3 (left panel). Reflecting the relatively young age 
a prominent, nearly vertical MS is present, together with 
the bright star at J w 6 and (J — H) w 0.55 (the red 
SG BD + 60°2532). Because of its low latitude, field stars 
(mostly disc) contaminate the CMD particularly at faint 
magnitudes and red colours. 

The intermediate-age nature of NGC 3960 is evident 
in the truncated MS and the presence of a clump of red 
evolved stars at (J — H) « 0.6 (right panel of Fig. |2J). 
However, conspicuous field-star contamination precludes 
further speculations. 

3.1. Differential reddening 

Previous work has shown that differential reddening is 
present in the fields of NGC 3960 and particularly M 52 
(Sect. Spatially variable reddening may introduce ad- 
ditional uncertainties in the analysis, particularly with re- 
spect to the intrinsic CMD morphology. To take this effect 
into account we divided the cluster fields into grids with 
cells of 5' width, along right ascension and declination, 
extending 17.5' in both directions. These dimensions cor- 
respond roughly to the cluster sizes (Sect. 13.41) . The 5' 
cell dimension provides both spatial resolution and sta- 
tistical significance in the number of stars in most cells. 
We illustrate this procedure in Fig. 01 for M52. For more 
accurate results we work only with the stars that remain 
after applying the colour-magnitude filter that discards 
stars with non-cluster colours (Sect. 13.4(1 . In this plot 
the projected stellar density of the central region of M 52 
clearly detaches from the outer regions. The correspond- 
ing J x (J — H) cell CMDs are shown in the right panel 
of Fig. □ 

Taking as reference the CMD of the central cell 
(AX = AY = 0') we built a fiducial line for the MS 
and evolved stars that represents the mean (J — H) colour 
distribution in bins of 0.25 mag in J. Because bright stars 
are scarce (particularly in M52), we used the respective 
Padova isochrone (Sect. 13.3(1 as a guide. The CMD fidu- 
cial lines of the central region of both clusters are shown 
in Fig. |21 For M 52 it is superimposed on all cell CMDs 
(right panel of Fig. 0) ■ To derive the differential reddening 
of a given cell with respect to the central one we applied 
a procedure that compares the observed CMD morphol- 

5 Following the Level 1 Requirement, according to 

http:/ /www.ipac.caltech. edu/2mass /releases /allsky /doc/ sec6=5al .html 
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Table 2. Differential reddening AE(B - V) in M52 
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-0.08 


+0.03 
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-0.05 
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-0.02 


0.01+0.07 


AE(B - V) 


0.07 + 0.08 


0.07 + 0.13 


0.07 + 0.09 


0.10 + 0.11 


0.11 + 0.12 


0.18 ± 0.14 


0.09 ± 0.07 





Table Notes. Uncertainties in AE(B - V) arc of the order of ±0.05. For absolute values: E(B - V) = 0.58 + AE(B - V). 
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Fig. 3. Left panel: spatial grid used to derive the differential reddening in M52. Only stars with colours consistent 
with the cluster's (Sect. 13 .4(1 are plotted here. North is up and East left. Right panel: J x ( J — H) CMDs extracted 
from the spatial grid (left panel). The fiducial line of the central cell is shown in all panels as a continuous line. Spatial 
orientation follows that in the left panel. AXo and AYo give the displacement along right ascension and declination 
of the center of each cell with respect to the central cell. 



ogy with the fiducial line. Differences are assumed to be 
caused only by reddening. The procedure varies E(B — V) 
from -0.8 to +0.8 in steps of 0.02. The observed J and H 
magnitudes of all stars in the cell are corrected for each 
new value of E(B — V), and the corresponding (J — H) 



colour is calculated. The chi-square is computed by sum- 
ming the squared difference of (J — H) colour between 
observed point and fiducial line over all the stars in the 
cell. The excess reddening (AE(B — V)) in each cell is the 
value of E(B — V) that produces the best match between 
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0.02 ± 0.09 
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Table Notes. Uncertainties in AE(B - V) are of the order of ±0.05. For absolute values: E(B - V) = 0.13 + AE(B - V). 



fiducial line and observed colours. The resulting values of 
AE(B - V) for M 52 are given in Table El where we in- 
clude as well the average values of AE{B — V) both for 
the right ascension and declination bins. 

As previously indicated (e.g. Danford & Thomas 1981) 
there is a reddening gradient from the south towards the 
north through the field of M52. With the present data 
we found that AE(B - V) ranges from -0.08 ± 0.05 to 
0.35 ± 0.05, and a least-squares fit to the data provides 
the relation AE(B - V) = (0.097 ± 0.030) + (0.0067 ± 
0.00319) x AY from south to north, with a correlation co- 
efficient CC = 0.98. With E(B -V) = 0.58 + 0.03 derived 
for the central region of M 52 (Sect. [£3, AE(B - V) val- 
ues convert to absolute reddening values of 0.50 ± 0.05 
and 0.93 ± 0.05, respectively. In right ascension, on the 
other hand, reddening is rather uniform, within uncer- 
tainties. The present values of the central and differential 
reddening agree with those of Danford & Thomas (JT981 ) , 
Kaltcheva l{T3971)> and Pandey et al. pTtTT)) . 

A similar procedure was applied to NGC 3960, and 
the respective values of AE(B — V) are given in Table [3] 
Contrarily to M52 the differential reddening is more 
uniform throughout the field of NGC 3960, within the 
uncertainties. AE{B — V) ranges from —0.10 ± 0.05 to 
0.21 ± 0.05 that, with the central value of E(B -V) = 
0.13 + 0.03 fSect. EO)) . convert to absolute reddening val- 
ues of 0.03 ± 0.05 and 0.34 ± 0.05, respectively. 

For the subsequent analyses we corrected the photom- 
etry for the differential reddening. We note, however, that 
the correction was applied only to the stars located at 
R < 20' from the cluster center. 

3.2. Field-star decontamination 

Field stars are usually an important component of wide- 
field CMDs, particularly of low-latitude star clusters 
and/or those projected against the bulge. Their presence 
in the fields of M 52 and NGC 3960 is clear (Fig.©, where 
they mimic an artificial sub-solar mass MS in M 52 and 



Table 4. Cluster and field-star statistics 



M52 NGC 3960 

AR N oos fci ffs N ons f c i ff B 

(') (stars) (%) (%) (stars) (%) (%) 

0) (2) (3) (4) (5) (6) (7) 

0-3 313 63 37 210 51 49 

3-6 601 43 57 434 28 72 

6-9 726 21 79 582 11 89 

9-12 874 8 92 803 9 91 

12-15 1053 2 98 927 100 



Total field 3567 20 80 2029 18 82 

Table Notes. Cols. 2 and 5: number of observed stars in the 
region. Cols. 3 and 6: fraction of member stars. Cols. 4 and 
7: fraction of field stars. Total field is R = 15' for M 52 and 
R = 12' for NGC 3960. 

a red sequence in NGC 3960. To better evaluate this ef- 
fect we show in the left panels of Fig. 01 the J x ( J — H) 
CMDs of the central 5' of M 52 (top panel) and NGC 3960 
(bottom). On both CMDs we superimpose the respective 
field-star contribution, extracted from the region 39.686 < 
R(') < 40 that corresponds to the same projected area. 
Most of the faint stars (J > 14.2) in M52 and the red 
ones ((J — H) > 0.6) in both clusters are probably field 
stars. 

To retrieve the intrinsic cluster-CMD morphology we 
use a field-star decontamination procedure previously ap- 
plied in the analysis of low-contrast (Bica & Bonatto 
2005), young embedded (Bonatto, Santos Jr. & Bica 
EIMjl . and young (Bonatto et al. I2006bjl OCs. In the 
present case we apply this procedure to the differential- 
reddening corrected photometry. As the offset field we take 
the region 20 < R(') < 40, that is large enough to pro- 
duce statistical representativity of the field-stars, both in 
magnitude and colours. The number of stars in both off- 
set fields is « 15 000. The algorithm divides the CMD 
in colour/magnitude cells from which stars are randomly 
subtracted in a fraction consistent with the expected num- 
ber of field stars. The remaining stars end up located in 
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(J-H) (J-H) 

Fig. 4. Panel (a): field stars (V) are superimposed on 
the differential-reddening corrected CMD of the central 5' 
of M52. Extraction areas match. Panel (b): field-star de- 
contaminated CMD. Solid line: 60Myr Padova isochrone 
with (m - M)j = 11.3 ±0.1 and E(J - H) = 0.18 ±0.01. 
Bottom panels: same as the top panels for NGC 3960. 
Solid line in panel (d): 1.1 Gyr Padova isochrone with 
(m - M)j = 11.3 ± 0.1 and E(J - H) = 0.04 ± 0.01. 
Colour-magnitude filters used to isolate MS/evolved stars 
are shown as dotted lines in panels (b) and (d). 

CMD regions where the stellar density presents a clear ex- 
cess over the field and consequently they have a high prob- 
ability of being cluster members. Because it actually ex- 
cludes stars from the original files - thus artificially chang- 
ing both the radial distribution of stars and LFs - we use 
field-star decontamination exclusively to derive the intrin- 
sic CMD morphology. Further details are given in Bonatto 
et al. l|2006bj) . 

In Table 0] we summarise the results of the decontam- 
ination procedure that quantify the presence of field stars 
throughout the fields of M52 and NGC 3960. Member 
stars appear to be found up to w 15' from the center 
of M52, while in NGC 3960 they extend to « 12'. The 
fraction of member stars in both fields is similar, ~ 20%. 
The field-star decontamination was applied to the whole 
range of colours and magnitudes (i.e. without discarding 
stars with non-cluster colours, Sect. 13. 3J) . Consequently, 



the above cluster dimensions are somewhat smaller than 
that suggested in Fig. for M 52 and those derived in 
Sect. 13.41 for both OCs undergoing colour-magnitude fil- 
tered photometry. This effect is better illustrated with the 
observed and colour-magnitude filtered radial density pro- 
files (Fig.®. 

The field-star decontaminated CMDs of the region 
R < 5' are shown in the right panels of Fig. where 
error bars are included to show the magnitude of the 
photometric uncertainties. As expected, most of the faint 
stars in M 52 and red ones in NGC 3960 were eliminated 
by the decontamination procedure. According to Fig. 
M52's CMD morphology is typical of a cluster younger 
than ~ 100 Myr with a well-developed MS spanning ~ 5 
magnitudes, while the old age of NGC 3960 is reflected in 
the low turnoff and giant branch. 

The colour spread in the MS of both clusters may be 
partly intrinsic and partly accounted for by residual red- 
dening, since we could not compute differential reddening 
for the offset field (R > 20'). However, differential redden- 
ing effects in the near-IR are smaller than in the optical, 
since E(J — H) « 0.3 x E(B ~V). To further minimise 
this effect we work with colour-magnitude filters broad 
enough to include most of the colour spread of the high- 
probability member stars (right panels of Fig.0J. 

3.3. Cluster age and distance from the Sun 

Cluster age is derived by means of solar- metallicity Padova 
isochrones (Girardi et al. 2002) computed with the 2MASS 
J, H and K s filters 6 . 

The red SG together with the rather well-defined 
low-MS of M52 reaching J « 14 constrain the age to 
60 ±10 Myr. Parameters derived from the isochrone fit are 
the observed distance modulus (to — M)j = 11.3 ± 0.1 
and colour excess E(J — H) = 0.18 ± 0.01, converting 
to E(B - V) = 0.58 ± 0.03. This age-solution is plotted 
in panel (b) of Fig. 0] The absolute distance modulus is 
(to - M) = 10.8±0.1, resulting in d & = 1.4±0.2 kpc. MS 
stars are restricted to the mass range 1.3 < m(M Q ) < 6.3. 
The Galactocentric distance of M 52 is dec = 8.7±0.2kpc, 
using Ro — 8.0 kpc as the Sun's distance to the Galactic 
center (Reid I1993J1 . However, with the recently derived 
value of Ro = 7.2 kpc (based on updated parameters 
of globular clusters - Bica et al. |2006), M52 is located 
dec = 7.9 ± 0.2 kpc from the Galactic center. M 52 is lo- 
cated fts 0.7 kpc outside the Solar circle. 

For NGC 3960 we derive an age of 1.1 ± 0.1 Gyr, 
(m-M)j = 11.3 ±0.1 and E(J - H) = 0.04 ±0.01, cor- 
responding to E(B - V) = 0.13 ±0.03, d e = 1.7 ±0.2 kpc 
and doc = 7.5 ± 0.2 kpc (Ro = 8.0 kpc) or dec = 
6.7±0.2kpc {R = 7.2kpc). Either way, NGC3960 is 
« 0.5 kpc inside the Solar circle. Because of the 2MASS 
faint- magnitude limit (Sect. 03), MS stars are detected in 



http:/ /pleiadi.pd. astro.it/isoc._photsys . 1/isoc.photsys. 01. html 2MASS 
transmission filters produced isochrones very similar to the Johnson ones, with 
differences of at most 0.01 in (J - H) (Bonatto, Bica & GiraTdi \!ffl%\ . 
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Table 5. Structural parameters. 
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Cluster 
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(stars pc - 2 


(stars pc -2 ) (p c ) 


(PC) 


(PC) 


(i) 


(2) 


(3) (4) 


(5) 


(6) 






Before differential-reddening 


correction 




M52 


3.2 ± 0.2 


49 ± 13 0.80 ±0.14 


8.0 ± 1.0 




NGC 3960 


4.4 ±0.1 


35 ± 9 0.69 ±0.12 


6.0 ± 0.3 








After differential-reddening 


correction 




M52 


3.7 ± 0.2 


45 ± 10 0.91 ±0.14 


8.0 ± 1.0 


13.1 ± 2.2 


NGC 3960 


4.8 ± 0.2 


36 ± 9 0.62 ±0.11 


6.0 ± 0.8 


10.7 ± 3.7 



Table Notes. We express King's profile as tr(-R) — a bg + CT 0K / (1± (R/Kcorc ) 2 ) ■ 
(Tqk and Hcore were allowed to vary while <7f, g was kept fixed. 



the range 0.95 < m(Mo) < 1.98. This age solution is 
shown in panel (d) of Fig. 0] 

3.4. Cluster structure 

Cluster structure was inferred with the radial density pro- 
file (RDP), defined as the projected number-density of 
MS/evolved stars around the center. The stars were se- 
lected by applying the respective colour-magnitude fil- 
ters (right panels of Fig. 0} to the differential-reddening 
corrected photometry (Sect. 13.1(1 . The use of colour- 
magnitude filters to discard foreground/background field 
objects was previously applied in the analysis of the OCs 
M67 (Bonatto & Bica . NGC 3680 (Bonatto, Bica 
& Pavani OrMjl . NGC 188 (Bonatto, Bica & Santos Jr. 
IHU), NGC 6611 (Bonatt o, Santo s Jr. & Bica l2UUfJJ> . and 
NGC 4755 (Bonatto et al. l2006bjl . To avoid oversampling 
near the center and undersampling for large radii, the 
RDPs were built by counting stars in concentric rings 
with radius AR = 0.5' for < R(') < 5, AR = V 
for 5 < R(') < 14, AR = 2' for 14 < R{') < 30 and 
AR = 4' for R > 30'. The residual background level 
corresponds to the average number of stars in the region 
20 < R(') < 40, resulting in a bg = 0.656 ±0.035 stars (')" 2 
and a bg = 1.204 ± 0.036 stars (')~ 2 , respectively for M52 
and NGC 3960. 

Fig.0shows the MS/evolved stars' RDP of both clus- 
ters. For absolute comparison between clusters the radius 
scale was converted to parsecs and the number-density of 
stars to stars pc -2 using the distances derived in Sect. 13.31 
The statistical significance of the RDPs is indicated by 
the relatively small lcr Poisson error bars. We also show 
in Fig. the RDPs produced with photometry prior to 
the colour-magnitude filtering. Clearly, the filtered pro- 
files present less fluctuation and reveal more of the cluster 
structure than the observed ones. In particular, the ob- 
served profiles tend to underestimate the cluster's exten- 
sion. 

Structural parameters were derived by fitting the 
colour-magnitude filtered RDPs with the two-parameter 
King (1966a) surface density profile, which describes the 
intermediate and central regions of normal globular clus- 
ters (King[l966bl Trager, King & Diorgovski UMS) . The 
fits were performed using a nonlinear least-squares fit rou- 
tine that uses errors as weights. To minimise degrees of 
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Fig. 5. Filled circles: differential-reddening corrected + 
colour-magnitude filtered radial density profiles (RDPs) 
of M52 (left panel) and NGC 3960 (right panel). Dashed 
lines: two-parameter King profile where the core size is in- 
dicated. Shaded regions: stellar background level. For com- 
parison purposes the RDPs prior to the colour-magnitude 
filtering are shown as empty circles. 



freedom in the fit the background level {(Jbg) was kept 
constant, corresponding to the residual values (see above). 
Parameters derived from the fit are King's central density 
of stars (oqk) and core radius (i? C ore)- The resulting pa- 
rameters are given in Table [S] and the best-fit solutions 
are superimposed on the colour-magnitude filtered RDPs 
(Fig.|SJ|. The limiting radius (R\\ m ) of a cluster can be es- 
timated by considering the fluctuations of the RDPs with 
respect to the background. i?i; m describes where the RDP 
merges into the background and for practical purposes 
most of the cluster stars are contained within Ry im . M52 
has core and limiting radii ~ 50% and ~ 33% larger than 
those of NGC 3960, respectively (Table EJ. 

We also give in Table parameters derived from the 
RDPs built with differential reddening uncorrected pho- 
tometry. The uncorrected values of uqk, Rcorc, and Rn m 
agree with the corrected ones, within uncertainties. 

The differential reddening corrected RDPs allowed us 
to compute the tidal radii of M 52 and NGC 3960 us- 
ing the three-parameter King (1962) model. The results 
are Etidal = 13.1 ± 2.2 pc (~ 1.6 x i?i im ) and iJtidai = 
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10.7±3.7pc (~ 1.8 x i?ii m ). The present ratios i?tidai/-Riim 
agree with those derived for bright OCs in Bonatto & Bica 
( 2005, that are in the range 1.4 — 1.9. 

The colour-magnitude filtering provided parameters of 
M 52 to be compared with those of Nilakshi et al. ( 2002 ) 
and Kharchenko et al. Ip005). The latter works do not 
apply such a filtering method. Probably because of the 
deep RDPs produced by the present method, our values 
of i?coro and i?i; m for M52 are about 50% smaller and 
larger, respectively, than those of Nilakshi et al. (2002) 
and Kharchenko et al. I|2UU5|I . 

King's profile provides an excellent analytical represen- 
tation of the stellar RDP of M 52 from the external parts 
to the core. Since it follows from an isothermal (virialized) 
sphere, the close similarity of the stellar RDP with a King 
profile may suggest that the internal structure of M52 
(particularly the core) has reached some level of energy 
equipartition after ~ 60Myr fSect.lO)). For NGC 3960, 
on the other hand, King's profile fails to reproduce the 
innermost region, particularly in the core where there is 
a marked excess density of stars (Fig. EJ). In this region 
the RDP increases almost as a straight line towards the 
center. This fact may suggest post-core collapse in this 
~ 1 Gyr OC (Sect. 13.3(1 . as seen in part of the globular 
clusters (Trager, King & Djorgovski 1995J). We will return 
to this point in Sect.EJ 

4. Luminosity and mass functions 

Luminosity and mass functions (</>(m) = g^J for the core, 
halo and overall cluster are derived for M 52 and NGC 3960 
following the methods presented in Bonatto & Bica (2005 
and references therein). To maximize the statistical signif- 
icance of field-star counts we take as the offset field the 
region 20 < R(') < 40, that lies > 5' and > 8' beyond the 
limiting radii of M52 and NGC 3960 (Sect. Eg) , 

To isolate MS stars we apply the colour-magnitude fil- 
ters (Fig. 0} to the differential-reddening corrected pho- 
tometries (Sect. 13.10 . resulting in the MS ranges 8.8 < 
J < 14.2 and 11.8 < J < 15.4, respectively for M52 and 
NGC 3960. The faint-magnitude limit of the MS in both 
clusters is brighter than that of the 99.9% Completeness 
Limit (Sect. © . 

To take the residual field-star contamination into ac- 
count we build LFs for each cluster region and offset field 
separately. J, H and K s LFs are built by counting stars 
in magnitude bins from the respective faint magnitude 
limit to the turn-off, for cluster and offset field regions. 
Magnitude bins are wider in the upper MS than in the 
lower MS to avoid undersampling near the turn-off and 
oversampling at the faint limit. Corrections are made for 
different solid angles between offset field and cluster re- 
gions. Intrinsic LFs are obtained by subtracting the offset- 
field LFs (Fig. U); they are transformed into MFs us- 
ing the mass-luminosity relations obtained from the re- 
spective Padova isochrones and observed distance modulii 
(Sect. 13- These procedures are applied independently to 
the three 2MASS bands, and the final MFs combine the J, 



E 



: v Wi M52 

r 

86 

6t 
<&- 

1 

- X = 1 48+0.09 


(a 


I 


X = 


NGC3960 

1 .07+0.22 


(d) : 

% 


: go 

X = 1.65+0.12 


(b 


) 


' f i 

X 


= 1 


1 J 

.26+0.26 


(e) : 

1: 


i 1 1 t 

X = 0.89+0.12 


1— 

<<=) 


-I — 

t 


— h 

ili 

x = 


w 

-0 


(f) : 

74+0.35 



1.3 2 3 4 5 6 1 1.21.3 1.5 1.7 2 

m (MJ m (M ) 



Fig. 6. Overall (top panels), halo (middle) and core (bot- 
tom) mass functions of M 52 (left panels) and NGC 3960 
(right). Dashed lines: fits with <f>(m) oc m~( 1+x \ 



H and K s MFs. MS mass ranges are 1.3 < m(M ) < 6.1 
and 0.95 < m(M Q ) < 1.98, respectively for M52 and 
NGC 3960. We note that the 1.3 Mq low-mass limit in 
M 52 is consistent under contraction arguments with an 
age of w 60 Myr. Although expected to populate the MS of 
a w 1 Gyr-old OC, stars less massive than 0.95 Mq cannot 
be detected in NGC 3960 with 2MASS mostly because of 
the increasing fraction of field-star contamination in that 
mass range (as already indicated by panel (d) of Fig.rjJ. 

Fig. shows the overall, halo and core MFs of both 
clusters together with the fits with the function <fi(m) oc 
m~( 1+x ', and the MF slopes are given in Tables and 

Both the overall and halo MFs of M52 with slopes 
X = 1.48 ± 0.09 and \ = L65 ± °- 12 are slightly steeper 
than Salpeter's l|1955|l IMF {% = 1.35). The core, however, 
is flatter with \ — 0-89 ± 0.12. Within uncertainties the 
present core, halo and overall MF slopes agree with those 
of Pandey et al. ( 2001 ) for the equivalent mass range. 

The overall MF of NGC 3960 with slope X = 1.07±0.22 
is flatter than Salpeter's, while that of the halo (x = 1.26± 
0.26) is similar. The core is very flat with \ = — 0.74±0.35. 
The present overall MF slope is significantly flatter than 
that of Prisinzano et al. (2004). In M52 and particularly 
NGC 3960, MF slopes increase from core to halo, which 
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suggests mass segregation in the core. We will discuss this 
point further in Sect. [5] 

4.1. Cluster mass 

Parameters derived from the core, halo and overall MFs 
of M 52 and NGC 3960 are given in Tables El and □ The 
number of evolved stars (col. 2) in each cluster region was 
obtained by integration of the respective field-star sub- 
tracted LF for stars earlier than the turnoff. Multiplying 
this number by the mass at the turn-off (m ~ 6 Mq 
in M52 and m « 2 Mq in NGC 3960) yields an esti- 
mate of the mass stored in evolved stars (col. 3). The 
observed number of MS stars and corresponding mass 
(cols. 5 and 6, respectively) were derived by integrat- 
ing the MFs in the mass ranges 1.3 < m(Af©) < 6.1 
and 0.95 < to(M ) < 1.98, respectively for M52 and 
NGC 3960. 

The total number (col. 7) and mass (col. 8) of the 
observed stars in M52 is computed by adding the cor- 
responding values of the evolved and MS stars. We also 
give in Table E| the projected mass density (col. 9) and 
volume mass density (col. 10). This total mass estimate 
for M52 should be taken as a lower limit, since we are 
not computing the low-mass stars still in the process of 
reaching the MS. On the other hand, upper-limits can 
be computed if we assume that stars less massive than 
1.3 Mq, although not detected in the 2MASS photometry, 
are present in the cluster in the form of PMS stars. We 
do this by extrapolating the MFs down to the H-burning 
mass limit, 0.08 Mq assuming Kroupa's l|2001J) universal 
Initial Mass Function (IMF), in which x — 0.3 ± 0.5 for 
the range 0.08 < m(M Q ) < 0.5. For masses in the range 
0.5 < m(Mg) < 1.3 we use our values of x which, within 
uncertainties are similar to the slope in Kroupa (2001) 
for the equivalent mass range. The extrapolated values of 
number, mass, projected and volume densities are given in 
cols. 7 to 10 of the third set of entries in TableEJ According 
to this extrapolation the core may contain about twice the 
observed mass, while the total cluster mass may be a frac- 
tion of ~ 3 larger than the observed value. 

Because of the 1.1 Gyr age, we estimate the total mass 
locked up in stars in NGC 3960 by taking into account 
all stars from the turnoff down to the H-burning mass 
limit, 0.08 Mq. For the core we extrapolate the flat MF 
down to 0.08 Mq. For the halo and overall cluster, on the 
other hand, we assume the universal Initial Mass Function 
(IMF) of Kroupa (|2UtJT)) . The results are given in Table □ 

We also give in Tables El and [7| parameters computed 
without the differential-reddening correction. The results 
are not very sensitive to this correction, at least in the 
near-IR. 

The present mass determinations for M52 are a fac- 
tor k 2 (observed) and « 6 (extrapolated) larger than 
that of Bruch & Sanders l|1983f) . The difference in the 
mass estimates can be accounted for by different analyti- 
cal methods. In the present case we have based our results 



on colour-magnitude filtered photometry, while Bruch & 
Sanders (1983) estimated mass according to an empirical 
relation that converts relative to absolute masses. 

5. Dynamical state 

5.1. M52 

Within uncertainties, the overall MF slope (x = 1-48 ± 
0.09) in the mass range 1.3 < m(M Q ) < 6.1 agrees with 
that of a Salpeter (x = 1.35) IMF. However, the MF slopes 
present significant spatial variations throughout the clus- 
ter, being flat (x = 0.89 ± 0.12) in the core and rather 
steep (x = 1-65 ± 0.12) in the halo (Table Eland Fig. Ell- 
in older clusters such spatial variation of MF slope can be 
accounted for by dynamical mass segregation, in the sense 
that low-mass stars originally in the core are transferred 
to the cluster's outskirts, while massive stars sink in the 
core. This process produces a flat core MF and a steep 
one in the halo (e.g. Bonatto & Bica 2005 and references 
therein) . 

Mass segregation in a star cluster scales with the re- 
laxation time t Te \ = 8 fa N t CT , where t cr = R/a v is the 
crossing time, N is the (total) number of stars and a v is 
the velocity dispersion (Binney & Tremaine[1987} . The re- 
laxation time is the characteristic time-scale for a cluster 
to reach some level of energy equipartition. For a typi- 
cal a v « Skins" 1 (Binney & Merrificld 1998) we obtain 
i rc i(overall) = 27 ± 4Myr and i rel (core) = 0.8 ± 0.2 Myr 
for the observed values; using the extrapolated values 
we obtain i re i (overall) = 311 ± 204 Myr and t re i(core) = 
4.5 ± 2.6 Myr. The « 60 Myr age of M52 (Sect.© corre- 
sponds to ~ 78 x < re i(core) (observed) or ~ 13 x i re i(core) 
(extrapolated). Thus, some degree of MF slope flattening 
in the core is to be expected. However, the ratio cluster 
age to t le \ drops to ~ 2.2 (observed) or ~ 0.2 (extrapo- 
lated) for the overall cluster, consistent with the Salpeter 
slope and absence of large-scale mass segregation. 

5.2. NGC 3960 

The flat core MF of NGC 3960 (x = -0.74 ± 0.35) 
seems consistent with the relaxation time t rc i(core) = 
1.1 ± 0.7 Myr, which corresponds to ~ 10 -3 of the clus- 
ter age. For the overall cluster we obtain t le \ (overall) = 
100 ± 63 Myr, corresponding to ~ 0.1 of the cluster age, 
which might explain its sub-Salpeter (x — 1.07 ± 0.22) 
MF. Large-scale mass segregation through the body of 
NGC 3960 is expected to have occurred, as well as low- 
mass star evaporation. NGC 3960 is dynamically more ad- 
vanced than M 52 (Sect.®}. 

6. Diagnostic diagrams of open cluster parameters 

Bonatto & Bica ( 2005 |) derived a set of parameters related 
to the structure and dynamical evolution of OCs in differ- 
ent dynamical states. We analysed nearby OCs with ages 
in the range 70-7 000 Myr and masses in 400 - 5 300 M , 
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Table 6. Parameters derived from the MFs of M52 



Evolved Observed MS Evolved + Observed MS stars 

Region N* m \ N* m bs N* map r 

(stars) (M ) (10 2 stars) (1O 2 M ) (10 2 stars) (1O 2 M ) (M pc -2 ) (M® pc" 3 ) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 

Before differential-reddening correction 
Core 0.73 ±0.11 0.9 ±0.1 2.2 ± 0.4 0.9 ±0.1 2.2 ± 0.4 109 ± 17 102 ± 16 97 ± 26 

Halo 3±1 18 ±6 1.69 ±0.11 4.7 ±0.6 9.7 ±1.2 4.7 ± 0.6 9.8 ± 1.3 4.9 ± 0.6 0.46 ± 0.06 
Overall 3±1 18 ± 6 1.53 ± 0.09 5.3 ± 0.6 11.3 ± 1.3 5.4 ± 0.7 11.5 ± 1.3 5.7 ±0.6 0.5 ± 0.1 2.2 ±0.5 

After differential-reddening correction 
Core 0.89 ±0.12 1.0 ± 0.2 2.4 ± 0.4 1.0 ±0.1 2.3 ± 0.4 91 ± 15 75 ± 12 78 ± 20 

Halo 3±1 18 ±6 1.65 ±0.12 4.6 ± 0.6 10.0 ± 1.4 4.7 ± 0.6 10.1 ± 1.4 5.1 ±0.7 0.47 ± 0.07 
Overall 3±1 18 ± 6 1.48 ± 0.09 5.5 ± 0.5 11.9 ± 1.3 5.4 ± 0.6 12.1 ± 1.4 6.0 ± 0.7 0.6 ± 0.1 2.2 ±0.6 

Evolved ± Observed + Extrapolated MS stars 
Core 0.89 ±0.12 1.0 ±0.2 2.4 ± 0.4 8.3 ± 5.4 4.8 ± 1.1 184 ±41 152 ± 34 13 ± 8 

Halo 3±1 18 ±6 1.65 ±0.12 4.6 ± 0.6 10.0 ± 1.4 82 ± 60 34 ±11 17 ±6 1.6 ± 0.5 
Overall 3±1 18 ± 6 1.48 ± 0.09 5.5 ± 0.5 11.9 ± 1.3 88 ± 64 38 ± 12 19 ± 6 1.8 ±0.6 0.2 ±0.1 
Table Notes. MS mass range: 1.3 — 6.1 Mq. Col. 4: MF slope. Col. 11: dynamical-evolution parameter r = age/t rc i. 

Table 7. Parameters derived from the MFs of NGC 3960 



Evolved 



Observed MS 



Evolved ± Observed + Extrapolated MS stars 



Region N* m 
(stars) (M Q ) 
(1) (2) (3) 



(4) 



N* TOobs N* map 

(10 2 stars) (10 2 M Q ) (10 2 stars) (10 2 M Q ) (M pc" 2 ) (M pc" 3 ) 
(5) (6) (7) (8) (9) (10) 



(11) 



Before differential-reddening correction 



Core 5±1 10 ± 2 -0.75 ± 0.31 0.3 ±0.1 0.5 ± 0.1 0.8 ±0.1 0.8 ±0.1 48 ± 6 
Halo 22 ± 5 44 ± 10 1.40 ±0.26 2.5 ± 0.3 3.3 ± 0.4 41 ±31 14.4 ±5.7 13 ± 5 
Overall 25 ± 6 50 ± 12 1.16 ±0.23 2.8 ± 0.3 3.7 ± 0.4 38 ± 27 14.1 ±5.1 12 ± 4 



52 ± 6 2300 ± 520 
1.6 ±0.6 

1.6 ±0.6 10 ±6 



After differential-reddening correction 



Core 5±1 10 ± 2 -0.74 ±0.35 0.3 ± 0.1 0.5 ± 0.1 2.5 ±1.7 1.1 ±0.3 94 ± 28 114 ±33 1010 ±604 
Halo 22 ± 5 44 ± 10 1.26 ± 0.26 2.3 ± 0.3 3.1 ±0.4 35 ± 25 12.6 ±4.7 11.2 ±4.2 1.4 ± 0.5 
Overall 25 ± 6 50 ± 12 1.07 ± 0.22 2.7 ± 0.3 3.6 ± 0.4 34 ± 24 13.0 ±4.5 11.5 ±4.0 1.4 ±0.5 11 ±7 



Table Notes. Same as Table El for NGC 3960. Detected MS mass range: 0.95 



;m< 



©■ 



following most of the present methodology. The evolution- 
ary parameter r = age/i re i (col. 11 of Tables and 01 
was found to be a good indicator of dynamical states. In 
particular, significant flattening in core and overall MFs 
due to dynamical effects such as mass segregation is ex- 
pected to occur for r corc > 100 and T ovcra n > 7, respec- 
tively. To these clusters were added the very young OC 
(age ~ 1.3 Myr) NGC 6611 (Bonatto, Santos Jr. & Bica 
2005)) and th e young (age ~ 14 Myr) NGC 4755 (Bonatto 
et al.lMhjl. 



6.1. Structural parameters 

Structural parameters are dealt with in Fig. Because 
the observed values of mass, density and evolutionary pa- 
rameters of M 52 are lower limits, we plot as well the cor- 
responding extrapolated values. For NGC 3960 we only 
consider the extrapolated values. The core density of 
NGC 3960, and particularly the extrapolated one of M 52, 
follow the trend presented by massive OCs for the cor- 
responding ages (panel a). A similar trend occurs for the 
overall density (panel b) , particularly for the extrapolated 
value for M52. Both clusters follow the massive OC re- 
lation involving overall mass and core radius (panel c). 
Panels (d) and (e) appear to suggest a correlation of i? core 
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Fig. 7. Relations involving structural parameters of OCs. 
Open circles: clusters more massive than lOOOAf©. 
Squares with a dot: m < IOOOMq. Dotted lines: least- 
squares fits to nearby clusters. Square: M52. Triangle: 
NGC 3960. Shaded areas: la borders of least-squares fits. 
Arrows indicate lower-limit estimates of mass, density and 
evolutionary parameter. Observed and extrapolated val- 
ues of mass, density and evolutionary parameter of M 52 
are shown. In panels (d) and (e) we show both determi- 
nations of (Iqc ■ 



and i?i; m with dcCi but the scatter of points precludes 
conclusions. M 52 and NGC 3960 fit in the observed cor- 
relation of core and limiting radii (panel f). About 2/3 of 
the sample OCs - including M 52 - seem to define a trend 
of increasing core radius with age (panel g), while the re- 
maining ones - including NGC 3960 - appear to follow a 
sequence of decreasing i? core with age. Similar trends oc- 
cur in the relation i?ii m with age (panel h). NGC 3960 fits 
in the relation of core and overall mass (panel i), while 
M 52 agrees at the la level. Further details on parameter 
correlations are given in Bonatto & Bica ( 2005 ). 



6.2. Dynamical parameters 

Dynamical-evolution parameters of M 52 and NGC 3960 
are compared with those of other OCs in different dy- 
namical states in Fig. |Hl The Salpeter-like overall MF 
slope of M 52 is consistent with the relation of Xoverall x 




10° 10' 10° 10° 10* 
Age (Myr) 



Fig. 8. Relations involving evolutionary parameters of 
open clusters. Symbols as in Fig. [7| Panels (f) and (i): 
fit of the linear-decay function \{t) = \ — t/tf . 



Af overa xi for massive clusters (panel a); the same is valid 
for NGC 3960 at the la level. The core MF slope of 
M52 is consistent with the expected one for its core 
mass (panel b) while the core of NGC 3960 seems signif- 
icantly flatter. Both clusters fit in the relations involv- 
ing Toverall X T corc (panel c), Xovorall X Xcore (panel d), 

Tcore x Xcorc (panel e), and r ovcral i x Xovcraii (panel h). The 
distribution of Galactocentric distances in terms of age 
(panel g) indicates a correlation, as suggested by Lynga 
( 1982 ). This correlation probably reflects the average tidal 
disruption effect on clusters, which increases for smaller 
Galactocentric distances (e.g. Boutloukos & Lamers 2003 
Lamers et al. [2005 ) . 

Panels (f) and (i) suggest a systematic flattening with 
age of the core and overall MFs, particularly for the less- 
massive OCs. The deviant object is the w 7 Gyr, massive 
OC NGC 188 whose orbit avoids the inner regions of the 
Galaxy for most of the time (Bonatto, Bica & Santos Jr. 
2005 J and references therein). The core and overall MFs 
of M 52 have slopes consistent with those of other mas- 
sive OCs of similar age. NGC 3960, on the other hand, 
has core and overall MFs significantly flatter than those 
expected for its mass, although they may be consistent 
with its age (see below). Except for NGC 188 the age de- 
pendence of the core and overall MF slopes can be pa- 
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rameterised by the linear-decay function x(<) = \ Q — t/tj, 
where Xo represents the MF slope in the early phases and 
tf is the flattening time scale. For the core MF we derive 
Xo = 0.83 ± 0.12 and t f = 0.9 ± 0.1 Gyr (correlation coef- 
ficient CC = 0.87); the overall values are Xo = 1-52 ± 0.08 
and t f = 2.4 ± 0.5 Gyr (CC = 0.79). These relations sug- 
gest that OCs in general are formed with flat core and 
Salpeter-like overall MFs. At cluster birth the core MF 
already is much flatter than the overall one, similarly to 
what was observed in e.g. NGC 6611 (Bonatto, Santos Jr. 
& Bica |2QQ6*)|. Such early core flattening may be partly 
linked to primordial processes related to molecular-cloud 
fragmentation. Because of the cumulative effects of mass 
segregation, additional core MF flattening becomes signif- 
icant after ~ 100 Myr and proceeds on an approximately 
linear time scale of ~ 0.9 Gyr (panel f). Flattening in the 
overall MF takes ~ 3 x longer to become evident and pro- 
ceeds at a slower pace (panel i) than in the core, which 
probably reflects the average time scale associated with 
low-mass evaporation, tidal stripping and encounters with 
giant molecular clouds. Compared to the overall cluster, 
harder physical conditions in the core such as high density 
and small relaxation time account for the fast evolution 
rate (Bonatto & Bica 1200'BI and references therein). 

The present MF flattening sequences - particularly the 
overall one - are defined mostly by the less-massive clus- 
ters, including NGC 3960. At least for ages < 3.2 Gyr, a 
linear-decay function provides a good analytical represen- 
tation of the flattening, both for the core and overall MFs. 
It remains to determine the behaviour of the flattening in 
older clusters, whether it keeps linear until dissolution into 
the field (particularly for the less-massive OCs) or satu- 
rates at some threshold value. A larger sample especially 
containing old OCs of any mass is required to (i) better 
constrain the analytical form of the less-massive OC flat- 
tening function, (ii) check whether the core and overall 
MFs of the massive OCs flatten in a similar way as those 
of the less- massive ones, and (Hi) compute the core and 
overall MF flattening time scales. 

7. Concluding remarks 

In this paper we analyzed the structure and distribution 
of stars in the OCs M 52 and NGC 3960 with J, H and K s 
2MASS photometry, restricted to stars with observational 
uncertainties £j,h,k s < 0.2 mag. To derive cluster param- 
eters we took into account the differential reddening and 
applied field-star decontamination and colour-magnitude 
filters. 

Compared to previous works we found differences in 
some parameters, such as cluster radius, mass and MF 
slopes. Most of those studies were based on spatial and 
magnitude-limited data. However, the discrepancies can 
be largely accounted for by the absence of a proper field- 
star decontamination. 

Reddening throughout the field of M52 varies from 
E(B - V) w 0.50 to E(B - V) w 0.93, with a marked gra- 
dient from south to north and a near uniformity in right 



ascension. For NGC 3960 we derive the range E(B — V) « 
0.03 — 0.34, uniform both in right ascension and decli- 
nation. The projected area considered in both fields was 
35' x 35'. 

For M52 we derived an age of 60 ± 10 Myr, absolute 
distance modulus (m — M)o = 10.8 ± 0.1, and distance 
from the Sun d & = 1.4 ± 0.2 kpc. For NGC 3960 we found 
an age of 1.1 ±0.1 Gyr, (to - M)o = 11.2 ±0.1, and d Q = 
1.7 ± 0.2 kpc. M52 is located « 0.7 kpc outside the Solar 
circle, while NGC 3960 is w 0.5 kpc inside it. 

Radial density profiles built with colour-magnitude fil- 
tered photometry are much deeper and yield more con- 
strained OC structural parameters than the observed 
ones. In particular, King's model provided an excellent fit 
to the differential-reddening corrected radial density pro- 
file of M52, from the external parts to the inner region. 
We found core and limiting radii of i? C ore = 0.91±0.14pc 
and R nm = 8.0±1.0pc. The inner region of NGC 3960, on 
the other hand, presents a clear excess of the stellar den- 
sity over King's profile, which may suggest post-core col- 
lapse in this 1.1 Gyr OC. Parameters derived are i? C orc = 
0.62 ± 0.11 pc and R hm = 6.0 ± 0.8 pc. The tidal radii 
of M 52 and NGC 3960, computed with the differential- 
reddening corrected RDPs are i?tidai = 13.1 ± 2.2 pc and 
i?tidai = 10.7 ± 3.7 pc. 

The present use of colour-magnitude filtered photome- 
try to build RDPs excludes a large fraction of non-cluster 
stars. As a result, observed and colour-magnitude filtered 
RDPs present significant differences (especially depth and 
shape) as shown in Fig. Because of the small resid- 
ual field-star contamination in the present RDPs, they 
can be used to derive more realistic intrinsic structural 
parameters. Nilakshi et al. lj2002H built RDPs using a 
photographic database (DSS photometry). They pointed 
out that the angular size was difficult to determine accu- 
rately due to the weak contrast between cluster and field 
stars. Kharchenko et al. Ij2005|) derived parameters based 
on stars with available proper motion, which restricts the 
analysis to brighter stars. For evolved clusters in which 
a significant fraction of low-mass stars has already been 
displaced from the center outwards, their brighter cutoff 
photometry may explain differences, together with lower 
statistics for a smaller sample of cluster stars. 

The mass stored in observed stars (1.3 < m(M ) < 
6.1) in M52 is TO ovcra n 1.2 x 1O 3 M ; the core contains 
Ti coro ~ 240 Mq. Extrapolating to the H-burning limit 
with Kroupa's IMF these values become m ovcra n 3.8 x 
10 3 Af Q and TO corc w 480 M Q . For NGC 3960 we derived 
^overall ~ 1-3 x 10 3 M Q and to coio w 114M Q . 

We conclude that both structurally and dynamically, 
the core and overall parameters of M 52 are consistent with 
those expected of an OC more massive than 1 000 M© 
and k 60 Myr old. The somewhat flat core MF slope 
(Xcoro ~ 0.89) as compared to the steeper one at the 
halo (xhaio ~ 1.65) and the Salpeter-like overall MF 
(Xovcraii ~ 1-48) suggests mild mass segregation in the 
inner region. This picture is consistent with the low val- 
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ues of the core/overall dynamical-evolution parameters, 
Tcorc ~ 13 - 78 and r overaU w 0.2 - 2.2. 

The parameters of NGC 3960, on the other hand, sug- 
gest that this OC - particularly the core - has reached an 
advanced dynamical state characterised by r corc ~ 10 3 
and T ovcra n ~ 11, respectively for the core and overall 
cluster. NGC 3960 presents evident signs of mass seg- 
regation in the core/halo region, with a very flat MF 
slope (xcorc ~ —0.74) and a steeper one at the halo 
(Xhaio ~ 1.26). Besides, the somewhat fiat (compared to 
Salpeter) overall MF (xovoraii ~ 1.07) suggests low-mass 
star evaporation through the cluster rim. Tidal disruption 
due to the Galactic field may have boosted the dynamical 
evolution of NGC 3960, since it is located w 0.5 kpc inside 
the Solar circle. 

Changes in parameters of M 52 and NGC 3960 because 
of differential-reddening correction are not significant, at 
least in the near-IR and/or because of the reddening 
ranges E(B - V) = 0.50 - 0.93 (M52) and E(B - V) = 
0.03 - 0.34 (NGC 3960). With respect to the analytical 
methods we conclude that colour-magnitude filters are es- 
sential to improve cluster parameters, particularly those 
related to structure. 
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